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ABSTRACT 



An improved dual tone mnltifrequency (DTMF) signal 
detector which uses the Goertzel DFT algorithm and which 
utilizes variable or differing frame widths that are frequency 
dependent for improved detection and reduced error. The 
DTMF detector includes a codec receiver which receives 
signals from the transmission media, and a digital signal 
processor (DSP) coupled to the codec. The DSP receives the 
digital samples and preferably applies the Goertzel DFT 
algorithm using differing frame lengths according to the 
present invention. The DTMF detector utilizes a different 
frame width for different tones of the possible tone frequen- 
cies according to the present invention. Thus the calculation 
uses a different frame length N for different ones of the 
uncorrected frequencies, wherein the different frame 
lengths comprise at least a subset N of the number of digital 
samples. The different frame lengths N are designed to 
optimally align the calculated frequency spectrum at each of 
the different uncorrelated frequencies. The calculation pro- 
duces an energy value for each of the different uncorrelated 
frequencies. The DTMF detector preferably multiplies a 
gain value with each of the energy values to adjust the gain 
of each of the energy values. After the frequency domain 
calculation, i.e., after energy values have been calculated for 
each of the different uncorrelated frequencies, and after any 
desired gain adjustment, the DSP determines maximum 
values of the energy values for each of the two or more 
frequency groups to detect the plurality of tones in the 
received signal. The DTMF detector also performs various 
other calculations to ensure valid tone detection. 

28 Claims, 10 Drawing Sheets 



ANALOG 
SIGNAL 
(DTMF+ 
SPEECH+ 
NOISE) 




DIGITAL 
SIGNAL 
SAMPLED 
AT 8KHz 
*(n) 



124 



GOERTZEL 
DFT USING 
DIFFERENT 
FRAME 
LENGTHS 



A(m) 
m«1,.16 



126 



GAIN 
ADJUSTMENT 



A'(m) 
m«1„16 



128 



FIND MAX 
GAIN IN 
EACH 
FREQUENCY 
GROUP 



M(m) 
m=1..4 


STATIC 
THRESHOLDING 


PASS 


DYNAMIC 
THRESHOLDING 
(SNR ESTIMATE) 


PASS 


SECOND 
HARMONIC 
THRESHOLDING 
(SPEECH 
DETECTION* 


PASS 


GUARD- 
TIME 
CHECK 














132 / J FAIL 




134 / 


FAIL 






FAIL 



,136 



I NO DETECTION I 



PASS 


TWIST 
COMPUTATION 
& 

THRESHOLDING 


PASS 


DECODE THE 
DIGIT USING 


DTMF 
DIGIT 






COLUMN & 
ROW INDEX 





FAIL 



FAIL 



138 



|N0 DETECTION 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul. 1, 1997 sheet 1 of 10 5,644,634 



DTMF 
DETECTOR 



102 



1 



103 




FIG. 1 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent jui. 1, 1997 sheet 2 of 10 5,644,634 




o 

1 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul i, 1997 sheet 3 of 10 5,644,634 



DIGITAL 
SIGNAL 
SAMPLED 
AT 8KHz 
x(n), n-1.~,178 



GDFT(1) 



GOERTZEL DFT 
N=N(1) 



GDFT(2) 



GOERTZEL DFT 
N=N(2) 



GDFT(3) 



GOERTZEL DFT 
N=N(3) 



GDFT(4) 



GOERTZEL DFT 
N»N(4) 



GDFT(16) 



GOERTZEL DFT 
N=N(16) " 



AO) 



A(2) 



A(3) 



A(4) 



A(16) 





g(D 

A'(D 



9(2) 
A'(2) 




9(3) 
A'(3) 




9(4) 
A'(4) 




g(ie) 

A'(16) 



N(1)=172; N(2)=177; N(3>=178; N(4)=178; N(5)=172 
N(6)=168; N(7)=168; N(8)=176; N(m)=N(m-8), for m=9...,l6 

FIG, 3 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent jui. 1, 1997 sheet 4 of 10 5,644,634 





Y k (N) - X(k) 




2cos(2*Pf*k/N) 



K 

"W N 



NOTE: 

FEEDFORWARD ONLY COMPUTED 
WHEN n=N FOR FINAL RESULT 

K 

-W N . - -exp( -J * 2 Pi * k/N) 
A(k) - amp X(k) 



FIG. 4 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul i, 1997 sheet 5 of 10 5,644,634 

151 





M(1)=0 
for0-1 rf 4) { 
if A'(i)>M(1) then 

endif 

} 


1(1). M(1) ., 


AY2) 


A'(3) 


1(2}. U(2) 


A'U) 


A'(5) . 


152 

/ 


M(2)«0 
for(i=5..8) J 
if AYi)>M(2) then 

l(2H 
endif 

} 


AYS} 


A*(7) 




AYS^ 


A*(9) t 


153 

/ 


M(3)=0 
for(i=9..12) { 
if A'(i)>M(3> then 

endif 

\ 


AY101 


A'(11) 




AY12} 


ATI 3). 


154 

/ 


M(4)=0 
for(l=13..16) { 
if A'(i)>M(4) then 

M(4)=A'(0 
endif 

\ 


A'(14) 


A 4 (15) 




A'(16) 



FIG. 5 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul i, 1997 



Sheet 6 of 10 



5,644,634 



Ts: STATIC THRESHOLD 



242 

M(1),M(2) YES 
^M(1)>Ts> Tt:5 



2U 



tl(2)>Ts 



YES 



CONTINUE TO DYNAMIC 



NO 


NO 


THRESHOLDING 






NO DETECTION 





FIG. 6 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul. 1, 1997 Sheet 7 of 10 



5,644,634 




12/11/2003, EAST Version: 1.4.1 



U.S. Patent 



JuL 1, 1997 



Sheet 8 of 10 



5,644,634 




12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul i, 1997 sheet 9 of 10 5,644,634 




YES 



Note: M'(0 IS THE PREVIOUS VALUE OF M! 
1(0 IS THE PREVIOUS VALUE OF 1(1). 
Tg: GUARD_TIME RATIO 



FIG, 9 



12/11/2003, EAST Version: 1.4.1 



U.S. Patent Jul. 1, 1997 Sheet 10 of 10 



5,644,634 




12/11/2003, EAST Version: 1.4.1 



5,644,634 

1 2 

SYSTEM AND METHOD FOR DUAL TONE allowance less than 1.5%. The term 'twist" refers to the 

MULTIFREQUENCY DETECTION USING difference, in decibels, between the amplitude of the stron- 

VARIABLE FRAME WIDTHS gest key pad row tone and the amplitude of the strongest key 

pad column tone. In general, the twist is required to be 

5 between -4 and +8 decibels. The noise immunity criteria 

FIELD OF THE INVENTION requires that if the signal has a signal to noise ratio (SNR) 

The present invention relates to the detection of dual tone n decibels, then the DTMF detector is required 

multifrequency coded signals, and more particularly to a t0 not ^ ^ i- e -» is required to detect the signal. 

DTMF detector which computes the frequency spectrum of ^ toe tf a tone a 

a received signal using a variable frame size for increased 10 duration greater than 40 inilliseconds, the DTMF detector is 

speed and accuracy. required to detect the tone, whereas if the tone has a duration 

less than 20 inilliseconds, the DTMF detector is required to 

DESCRIPTION OF THE RELATED ART not detect the tone. Speech immunity refers to the ability of 

Dual tone multifrequency (DTMF) coding is a generic * e ^ ***** t0 <H*bg** ^TMF tone 

name for push-button telephone signaling which is used in 15 T^^^t 8 ^ 

North American telephone systems A DTMF signal is used ^ * abdl * of ^ ™* to Y™* 

for transmitting a phone number or the like from a push ff** s m t &P £ 'J^> ^T^ 0 ? 

button telephone to a telephone central office. DTMF sig- ^/™^™ e ^ ^ °T ^xt^ 0 ™? 

naling is quickly replacing dial pulse signaling in telephone ^ertently detectmg speech as a DTMF signal, 

networks worldwide. In addition to telephone call signaling, 20 Originally, DTMF signal generators and detectors utilized 

DTMF coding is also becoming popular in interactive con- analog circuitry to generate and decode DTMF signals, 

trol applications, such as telephone banking, electronic mail However, with the rapid advance of VLSI technology and 

systems, and answering machines, wherein the user can digital signal processing (DSP) technology, many DTMF 

select options from a menu by sending DTMF signals from systems are now employing digital signal processors for 

a telephone. 25 increased accuracy and cost efficiency. The advantages of a 

Adi^tonesignalisrer^esentedbytwosinusoidalsignals ^ ^ IMF generation and detection system include 

whosefrequendesareseparatedmbandwidthandwhichare ™P™ed accuracy, precision, stability, versatility and 

uncorrected in order to avoid false tone detections. In reprograjnmability, as well as lower chip^unt The advan- 

general, normal speech patterns or noise signals produce a 30 ^ es * a DSP implementation of a DTMF system are 

signal with energy distributed throughout the frequency 30 e T CU ^^ * telep 1 hon . c ^^Cented Office, 

band. In some cases, speech may have sufficient energy on where DTMF detection can be simultaneously performed on 

both DTMF frequencies to trigger a false detection. mum P le tele P hone channels. 

As noted above, a DTMF signal encoder generates a ^ g eneral i a DTMF detector examines the line or com- 
DTMF signal by adding together two sinusoidal signals. A 35 ™™ cation channel for the presence of two sinusoids using 
DTMF signal includes one of four tones, each having a dedicated frequency domain algorithms, including modified 
frequency in a low frequency band, and one of four tones, Goertzel algorithms, DFT/FFTs, auto-correlation, zero 
each having a frequency in a high frequency band. In current crossing counting, and narrow band filter-based methods, 
DTMF systems, the following frequencies are allocated for others - ^ Goertzel algorithm is the most popular 
the four tones in the low frequency band: FA=697 Hz; 40 ^thm used in DTMF signal detection because it has the 
FD=770 Hz; FC=852 Hz, and; FD=941 Hz; and the follow- accuracy of a DFT-based algorithm while also having com- 
ing frequencies are allocated for the four tones in high pactional efficiency comparable to narrowband filter-based 
frequency band: FE=1209 Hz; FF=1336 Hz; FG=1447 Hz, algorithms. Conventional use of the Goertzel algorithm is 
and FH-1633 Hz. The frequencies used for DTMF encoding similar t0 DFT methods, which compute values in the 
and detection are defined by the CdTT and are accepted 45 frequency domain at desired discrete points. As a result, this 
around the world, thus allowing dialing compatibility method has ^ disadvantage of the DFT, whereby the 
throughout the word. computed results may not be sufficiently close to the desired 

AsDmFsignalstmvddowna^^ 

a telephone line, the signals may become distorted due to ™°^ ThlS F°blem ^exacerbated by the fee mat the 

attenuation or to any number of other affects, such as 50 ™ f?™ d ^ tobt are deliberately 

channel noise, radiation, etc In addition, in some instances, chosen to * UDCOrrdated for detectlon - 

voice or speech signals are propagated down a telephone Therefore, in order to achieve sufficient resolution and 

line simultaneously with DTMF signals. Due to noise affects ^us sufficient accuracy, the window size or frame size is 

or speech, a DTMF signal can be missed by the detector. required to be fairly large. As a result, conventional DTMF 

This is referred to as "talk-down". Also, speech can have a 55 detectors using a fixed frame size to compute frequency 

frequency domain content similar to a DTMF signal and domain values generally also have the disadvantage of slow 

trigger erroneous detection, this being referred to as a detection and high computational costs to produce each 

"talk-off" effect. Noise other than speech can also cause Tesult lbas improved DTMF signal detection techniques 

adverse effects, and the receiver thus is sometimes unable to ^ c desired which have increased accuracy and improved 

discern which DTMF signal has been detected. 60 cmcienc y* 

„ ^"S^f ST TP? 4 V f° US ^ SUMMARY OF THE INVENTION 
dards for DTMF detectors. These standards involve various 

criteria, such as frequency distortion allowance, twist The present invention comprises an improved dual tone 

allowance, noise immunity, guard time, talk-down, talk-off, multifrequency (DTMF) or multitone signal detector which 

acceptable signal to noise ratio, and dynamic range, etc The 65 more efficiently and reliably detects DTMF signals. The 

distortion allowance criteria specifies that a DTMF detector present invention uses a frequency domain detection 

is required to detect a transmitted signal that has a distortion method, such as a modified version of the Goertzel DFT 
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algorithm, and utilizes variable or differing frame widths 
that are frequency dependent for improved detection and 
reduced error. 

The DTMF detector according to the present invention 
preferably includes a coder/decoder (codec) receiver which 
receives signals from the transmission media. The codec 
samples the received analog signals and produces digital 
signals, i.e., digital samples. The DTMF detector also 
includes a digital signal processor (DSP) coupled to the 
codec. The DSP receives the digital samples and preferably 
applies the Goertzel DFT algorithm using differing frame 
lengths according to the present invention. The present 
invention further includes a memory coupled to the DSP 
which is used by the DSP for temporary storage of data 
during processing (RAM), and retrieval of data such as filter 
coefficients (ROM). 

According to the present invention, the DTMF detector 
receives a plurality of digital samples of a received signal, 
wherein the received signal may include a plurality of tones, 
such as a DTMF or MTMF (multiple tone multifrequency) 
signal. The received signals may include DTMF or MTMF 
tone signals and also may include one or more speech 
signals and/or noise. The plurality of tones comprise two or 
more tones from a plurality of different uncorrelated fre- 
quencies. In one embodiment, the plurality of different 
uncorrelated frequencies comprise two or more frequency 
groups. After receiving the digital samples, the DTMF 
detector calculates a frequency spectrum of the plurality of 
digital samples for each of the plurality of different uncor- 
related frequencies. 

The DSP preferably uses the Goertzel algorithm to com- 
pute the frequency spectrum. The Goertzel DFT algorithm 
divides the sampled time domain signal into a plurality of 
discrete blocks or frames and then performs Fourier tech- 
niques to obtain energy values in the frequency domain. The 
DTMF detector utilizes a different frame width for different 
tones of the possible tone frequencies according to the 
present invention. Thus the calculation uses a different frame 
length N for different ones of the uncorrelated frequencies, 
wherein the different frame lengths comprise at least a subset 
N of the number of digital samples. The different frame 
lengths N are designed to optimally align the calculated 
frequency spectrum at each of the different uncorrelated 
frequencies. It is noted that the same frame length may be 
used for certain of the frequencies. The calculation produces 
an energy value for each of the different uncorrelated 
frequencies. 

It is noted that calculating the frequency spectrum using 
different frame lengths N produces differing energy contents 
in a plurality of the energy values. Thus the DTMF detector 
preferably multiplies a gain value with each of the energy 
values to adjust the gain of each of the energy values. After 
the frequency domain calculation, i.e., after energy values 
have been calculated for each of the different uncorrelated 
frequencies, and after any desired gain adjustment, the DSP 
determines maximum values of the energy values for each of 
the two or more frequency groups to detect the plurality of 
tones in the received signal. The DTMF detector also 
performs various other calculations to ensure valid tone 
detection. 

Therefore, the present invention comprises a DTMF 
detector and method for more accurately detecting the 
presence of two or more tones in a received signal. Hie 
present invention uses frequency domain techniques and 
uses different frame widths for different tones of the DTMF 
tone frequencies for improved speed and accuracy. 



BRIEF DESCRIPTION OF THE DRAWINGS 

A better understanding of the present invention can be 
obtained when the following detailed description of the 
5 preferred embodiment is considered in conjunction with the 
following drawings, in which: 

FIG. 1 is a block diagram of the DTMF detector of the 
present invention. 

FIG. 2 is a flowchart diagram illustrating operation of the 
10 DTMF detector of the present invention; 

FIG. 3 illustrates a plurality of filters in the DTMF 
detector of FIG. 1, which apply the Goertzel algorithm using 
a variable frame length according to the present invention; 
FIG. 4 illustrates operation of the Goertzel algorithm; 
FIG. 5 is a block diagram illustrating logic which finds the 
maximum gain in each group of detected tones; 

FIG. 6 illustrates operation of static thresholding in the 
flowchart diagram of FIG. 2; 
20 FIG. 7 illustrates operation of the dynamic thresholding in 
the flowchart diagram of FIG. 2; 

FIG. 8 illustrates operation of the second harmonic energy 
threshold check in the flowchart of FIG. 2; 
FIG. 9 illustrates operation of guard time checking func- 
25 tions performed by in the flowchart of FIG. 2; and 

FIG. 10 illustrates operation of twist checking functions 
in the flowchart of FIG. 2. 
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DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 



Incorporation by Reference 

The following U.S. patents and references are hereby 
incorporated by reference: 
35 U.S. Pat No. 5,408,529 titled ■'Dual Tone Detector Operable 
in the Presence of Speech or Background Noise and 
Method Therefor" — Greaves, issued Apr. 18, 1995, is 
hereby incorporated by reference. 
U.S. Pat. No. 5,428,680 titled "DTMF Signal Receiving 
40 Apparatus Equipped With a DTMF Signal Judging 
Circuit" — Murata et al., issued Jun. 27, 1995 is hereby 
incorporated by reference. 
U.S. Pat. No. 5,257,309 titled **Dual Tone Multifrequency 
Signal Detection and Identification Methods and 
45 Apparatus"— Brandman et al., issued Oct 26, 1993 is 
hereby incorporated by reference. 
U.S. Pat No. 5,353345 rifled "Method and Apparatus for 
DTMF Detection"— Galand, issued Oct 4, 1994 is hereby 
incorporated by reference, 
so U.S. Pat No. 5,426,696 titled "Method of Improving 
Receiver Sensitivity and Speech Immunity With DTMF- 
Reception" — Zimbrek, issued Jun, 20, 1995 is hereby 
incorporated by reference. 
U.S. Pat No. 5325,427 titled "Apparatus and Robust 
55 Method for Detecting Tones" — Dighe, issued Jun. 29, 
1994 is hereby incorporated by reference. 
"Digital Signal Processing Applications Using the ADSP- 
2100 Family", from Analog Devices Corporation, Volume 
1, Chapter 14, pages 441-500 is hereby incorporated by 
60 reference in its entirety. 
DTMF Detector 

The present invention comprises a DTMF detector which 
utilizes a variable frame size or differing frame sizes accord- 
ing to the present invention. The DTMF detector is com- 
65 prised in a digital telephone answering machine in the 
current embodiment However, it is noted that the present 
invention may be used in various other applications, includ- 



12/11/2003, EAST Version: 1.4.1 



5,644,634 

5 6 

ing telephone switching and interactive control applications, The DTMF detector 102 also preferably includes a digital 

telephone banking, fax on demand, etc. The DTMF detector signal processor (DSP) 106 coupled to the codec 104. The 

of the present invention may also be used as the DTMF DSP 106 is preferably from the ADSP-2100 family from 

detector in a telephone company Central Office, as desi red. Analog Devices, Inc. Other equivalent DSPs are acceptable, 

The present invention was designed to meet CCTIT and 5 although a DSP that can perform 16x16 bit hardware mul- 

AT&T standards, although the invention can be repro- tiplication is preferred for accuracy reasons. It is noted that 

grammed or undergo minor modifications to meet other a software multiplier may be used, although the operating 

telecommunication standards. The preferred embodiment d ^ ^ ^^^y reduced, 

was developed to perform DTMF detection functions and ^ DSp 106 recdves ^ les and preferably 

preferably compnses a dual tone detector which aetecte ftvo Hes fr dQm ^ techm ^ preferably me Goert- 

tones, wherein each tone is one of a plurality of predeter- tTYm* i -*u • • ui Ja *L % 

mined tones from two respective frequency groups. *7 algorithm, using variable ordiffenngframekngms 

However, the present invention may comprise a multi tone ^ccoroing to die preset invention. The Goertoel DFT algo- 

multi frequency (MTMF) detector for detecting a greater nthm divides the sampled time domain signal mto a plurahty 

number of different uncorrected tones in a received signal, of discrete blocks or frames and ^ performs Fourier 

Le., two or more different uncorrected tones. The two or 15 techniques on each frame to obtain energy values in the 

more different uncorrected tones should be spaced suffi- frequency domain for each of the possible tones. The DTMF 

ciently apart for reliable detection. The two or more tones detector 102 utilizes differing frame widths far different 

may be from two or more different frequency groups or from ones of the possible tone frequencies according to the 

a single frequency group, as desired. Where two or more present invention. In other words, the DTMF detector 102 

tones may be transmitted from a single frequency group, the 20 utilizes differing frame widths for at least a plurality of the 

possible frequencies or tones in this single frequency group possible tone frequencies. 

comprise different uncorrelated frequencies that are each The DTMF detector 102 further includes a memory 108 

spaced sufficiently apart such that any two or more tones coupled to the DSP which is used by the DSP 106 for storage 

from the single frequency group can be transmitted with and retrieval of data. Preferably, the memory comprises a 

reliable detection. Thus the system and method of the 25 RAM (random access memory) for processing storage and 

present invention may be used for detecting any number of loading of program at the time of operation, and a ROM 

tones in a received signal. (read-only memory) for storage of fixed values, such as filter 

FIG. 1 — Block Diagram coefficients and other parameters used by the program. Hie 

Referring now to FIG. 1, a block diagram of a DTMF program is also preferably stored in the ROM before loaded 

detector 102 according to the present invention is shown. As 30 into RAM. The memory requirement for the DTMF detector 

shown, the DTMF detector 102 preferably comprises a port alone is about 2.0 kbytes of ROM and 200 bytes of RAM. 

or connector or other means 103 for receiving analog or The memory requirement is architecturally dependent and 

digital signals. The port 103 is adapted for coupling to a also depends on implementation. It is noted that the memory 

communications medium, such as a phone line, cable or 108 may comprise any of various types or sizes, as desired, 

other transmission line. It is noted that the port 103 may 35 It is noted that the present invention may include a general 

comprise any of various means or connectors for coupling to purpose microprocessor instead of DSP 106. The present 

a communications line. invention may also include dedicated digital and/or analog 

The DTMF detector 102 also preferably comprises a hardware instead of, or in addition to, the DSP 106. Thus, 

coder/decoder (codec) receiver 104 coupled to the port although the following description describes the DSP 106 

which receives an analog signal and converts the analog 40 performing the present invention, it is noted that the present 

signal into digital format The codec receiver 104 preferably invention may be performed by a general purpose 

samples the analog signal at 8 kHz and utilizes pulse code microprocessor, or the present invention may be at least 

modulation (PCM) or other suitable techniques to produce partially or totally implemented in digital or analog logic, as 

corresponding digital data. The codec chip 104 preferably desired The use of a programmable DSP 106 is the preferred 

comprises linear analog to digital (A/D) converters and 45 embodiment of the invention, 

digital to analog (D/A) converters. The codec receiver 104 FIG. 2 — Flowchart Diagram 

preferably comprises all the necessary A/D, D/A, sampling Referring now to FIG. 2, a flowchart diagram illustrating 

and filtering circuitry for bi-directional analog digital inter- operation of the DTMF detector 102 according to the 

facing. Once analog to digital conversion has been preferred embodiment of the present invention is shown, 

performed, digital data or digital samples are generated 50 The DTMF detector 102 is coupled to a transmission media, 

based upon the analog signal. As shown, the DTMF detector 102 receives an analog signal 

It is noted that the received signal may be compressed or and determines whether dual tone multi frequency (DnTMF) 
companded, and thus digital data produced by the A/D signals are comprised within the received signal. As shown, 
converter in the codec 104 may be companded, Le., may the received analog signal may comprise one or more of a 
comprise logarithmically compressed digital data. As is 55 DTMF signal, a speech or voice signal, or noise. In step 122, 
well-known in the art, companding refers to logarithmically the codec 104 in the DTMF detector 102 receives the analog 
compressing a signal at the source and expanding the signal signal and converts the analog signal into digital format In 
at the destination to obtain a high end-to-end dynamic range step 122 the codec receiver 104 preferably samples the 
while reducing dynamic range requirements within the com- analog signal at $ kHz and performs pulse code modulation 
munication channel. In this instance, the codec receiver 104 60 (PCM) or other suitable techniques to produce correspond- 
logarithmically expands the data to a linear format, prefer- ing digital data. As shown, the codec receiver 104 produces 
ably a 16 bit linear format. a digital signal referred to as x(n), which is sampled at 8 

In one emrxxiiraent, the port 103 in the DTMF detector kHz. 

102 receives digital signals (i.e., linear or logarithmic PCM) It is also noted that the received signal may be logarith- 

directly from the communication channel or transmission 65 mically compressed or companded, and thus step 122 may 

media, and thus a codec 104 is not required in the DTMF further comprise logarithmically expanding the compressed 

detector 102 for analog to digital conversion. data to produce linear data. 
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The following comprises a brief description of steps 

124-140, and a more detailed description of each of the TABLE 1 
steps follows. In step 124 the DSP 106 receives the digital 
signal x(n) from the codec 104 and performs the Goertzel 
DFT using differing frame lengths according to the present 5 
invention. In step 124 the DSP 106 produces a plurality of 
energy values A(m) where m=l-16. The plurality of energy 
values A(m) comprise energy values at each of the two sets 
of four frequency tones for the first and second harmonics of 
the frequencies. In step 126 the DSP 106 adjusts the gain of 10 

each of the energy values A(m) to compensate for differing ^ J . _ ^ ^ , 
energy content due to the differing frame lengths. The DSP ^ detection <rf the dual tones is performed by math- 
106 in step 126 produces a plurality of adjusted energy transferring the input time domain signal into 
values A'(m). In step 128 the DSP 106 determines the the frequency domam usmg vanous Founer transform tech- 
maximum level of the energy values A'(m) in each fre- is ^".^ * w ^°™> «*> ^crete Fourier transform 
quency group (DFT) is commonly used to transform discrete time domain 

In the following steps 130-138 the DSP 106 performs mt0 ^ ^ scrttc domai * components. As 

various tests to ensure valid tone detection. In steps 130 and ? so note( J m background section, the Goertzel algorithm 

132 the DSP 106 performs static and dynamic thresholding " a Pppular method for performing the DFT computation 

respectively to ensure that the energy values A*(m) meet 20 ™ e ? oertzel flg°nthm offers the accuracy of a DFT-based 

certain basic criteria. The static and dynamic thresholding calculation, while having increased computational efficiency 

performed in steps 130 and 132 eliminate invalid DTMF 1S com P arable to narrow band ^tei-based algorithms, 

tones based on both signal level and signal/noise ratio. After . * V 0 ** conventional usage of the Goertzel algo- 

static and dynamic thresholding are performed in steps 130 nthm is siimlar t0 DFr method ** d 15 Performed by 

and 132, in step 134 the DSP 106 performs second harmonic 25 com P uti *g frequency domain values at the desired discrete 

thresholding to prevent detection triggered by speech. Hie P° mts ' As a rcsult > me Goertzel algonthm also has the 

DTMF detector of the present invention examines the sec- disadvantage of the DFT, whereby the computed results may 

ond harmonics of the fundamental DTMF frequencies using not be sufficientl y close to the desired frequencies desired to 

a novel ratio comparison method according to the present * deeded, thus severely reducing accuracy. This is referred 

invention, which further distinguishes speech signals from 30 to as me leaka S e effecL 111 other words > me D™* 7 detector 

DTMF signals. ^ 2 ^ desi 8 n e<i to detect dual tones which will occur at 

In step 136 the DSP performs a guard time check to s P ccified predefined frequencies. In the preferred 

evaluate the frequency domain results and to ensure that the embodiment, these frequencies are 697 Hz, 770 Hz, 852 Hz, 

signal lasts at least a certain amount of time. Hie guard time 941 ^ 1209 Hz > 1336 Hz, 1477 Hz, and 1633 Hz. Hie 

check evaluates the maximum energy values in relation to 35 frequencies used in the DTMF standard are not spaced at 

prior and subsequent frames for improved detection. In step e( * uai mcre ments but rather are deliberately designed to be 

138, the DSP performs twist computation and thresholding. ^correlated f°f more accurate detection. 

Hie present invention uses a novel twist computation tech- As a result, if the conventional method of afixedN(frame 

nique which only performs the twist computation when the ske > is used m ^ Goertzcl ^g orithm > » e frequencies are 

received signal is deemed stable, thus increasing accuracy. 40 not optimally aligned at the frequency bins. The effect is that 

After the twist computation is performed in step 138, in step ±e desured outouts m not exacti y 31 me frequency bins 

140 the DSP decodes the digit using the column and row (multiples of Fs/N). This causes the energy to be distributed 

indices and outputs the pressed DTMF digit from the amon 8 n eighbonng frequency bins, referred to as leakage, 

received DTMF signals. As noted in FIG. 2. if any of the severely degrade the accuracy of the frequency 

steps 130-140 fail, then no detection is indicated. 45 domain ou ^ uts * ^ foUowill g toble illustrates the effect of 

Step 124-Goertzel DFT with Variable Frame Lengths leaka S e usm £ 311 example of the Goertzel algorithm for 

Referring again to FIG. 2, in step 124 the DSP 106 011417 frequency 852 Hz (N=205). 
receives the digital signal x(n) from the codec 104 and 

performs the Goertzel DFT for each of the possible tone TABLE 2 

frequencies using varying or differing frame lengths accord- 50 The effect of leakage 

ing to the present invention. As discussed in the background — — 

section, decoding a DTMF signal involves detecting two error ' m K O-Qi 

tones in the received signal and then determining the number wo 

pressed by the user at the sending end based on the values g Jn<aB) 

of the detected tones. More generally, decoding a MTMF 55 output 

signal involves detecting a plurality of tones in the received error(dB) 



0.05 


0.1 


0.5 


IX) 


5.0 


80.6 


74.7 


60S 


49.6 


43.6 


>13.4 


>19.3 


>33.1 


>44.4 


>50.4 



CQTT recommendations Q.23 and Q.24 in the <( Red 0t6 ' 

Book**, Volume VI define two groups of frequencies that are In general, it has been found that a fixed-frame size of 205 

used for in-band signaling. Group 1 comprises the frequen- 60 points produces the least error in a Goertzel/DFT computa- 

cies 697 Hz, 770 Hz, 852 Hz, and 941 Hz, and these tion. The inherent property of the DFT requires FNT to be 

frequencies identify the rows of the telephone keypad, and equal to one, where F is the frequency domain resolution, N 

Group 2 comprises the frequencies 1209 Hz, 1336 Hz, 1477 is the frame size in number of samples, and is also the 

Hz and 1633 Hz, and these frequencies identify the columns computed number of outputs in the frequency domain, and 

of the telephone keypad. The DTMF signals each contain 65 T is the smallest period and equals 1/Fs, where Fs is the 

one row and one column frequency, and these signals are sampling frequency. In telephony applications, the period T 

assigned using the following table. is equal to 1/Fs= 1/8000. Urns, in order to achieve sufficient 
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resolution, the frame size N is required to be very large. As 
a result, detection of speech is generally slow and compu- 
tational costs are higher for each result Even worse, a large 
frame size may result in the detection speed not meeting 
telecommunications standards. 

Therefore, prior art DTMF detectors which utilize the 
Goertzel algorithm and/or the DFT and which utilize a fixed 
frame size N only compute frequency domain values in 
multiples of Fs/N, wherein Fs is the sampling frequency, 
typically 8000 Hz for telephony applications. Therefore, for 
example, if Fs=8000 Hz and N=100 sample points, the 
Goertzel algorithm only computes results at 80 Hz, 160 Hz, 
240 Hz, etc. Therefore, in summary, since the DTMF 
standard frequency values are designed to be uncorrelated, 
in general, no single frame size value N can be chosen which 
' results in all of the desired frequencies lined up at the bins. 

Therefore, the present invention performs the Goertzel 
algorithm utilizing a varying or differing frame size N for a 
plurality of the DTMF frequencies to provide increased 



algorithm performs a second order recursive computation of 
the DFT by computing a new Yk(n) output for every new 
input sample x(n). The DFT result, referred to as X(k) is 
equivalent to Yk(n) when n=N, i.e., X(k>=Yk(N). As shown 

s in FIG. 4, the operation of the Goertzel algorithm can be 
divided into two phases. The first phase involves computing 
the feedback legs in FIG. 4. The second phase computes the 
feed forward when n=N and thus evaluates the value X(k). 
As discussed above, in the present invention the DSP 106 
performs the Goertzel algorithm shown in FIG. 4 using a 

10 different frame length for different ones of the DTMF 
frequency values. 

After the Goertzel DFT is applied to the received digital 
signal using variable frame lengths according to the present 
invention in step 124, the result is a plurality of frequency 
spectra values A(m) corresponding to the 8 DTMF standard 
frequencies or tones for the first and second harmonics. 

The following two tables compare the Goertzel algorithm 
implementation using a fixed frame size and a variable 
frame size. The first table below illustrates an example using 
N=205 as the fixed frame size and shows the K errors. It is 



accuracy and efficiency. The DTMF detector 102 of the 20 noted that a fixed frame size of N=205 is the best possible 



present invention preferably uses a different frame size for 
different DTMF tone frequencies and then adjusts the gain 
after the frequency domain values have been computed. This 
provides increased accuracy and speed over prior art meth- 
ods. The DTMF detector 102 of the present invention 
provides better frequency distortion allowance and better 
twist allowance. 

Referring again to FIG. 2, in step 124 the DSP 106 applies 
the Goertzel DFT method using different frame lengths 
according to the present invention. As shown, the DSP 106 
receives a digital signal referred to as x(n), wherein the 
digital signal x(n) comprises a plurality of samples of a 
received signal Referring now to FIG. 3, The DSP effec- 
tively computes 16 Goertzel DFTs with different frame 
lengths N(l)-N(16). The different frame lengths for each of 35 
the Goertzel DFTs are as follows: 



selection and produces the least maximum error for all K 
values. 



25 



30 



K = N » f/fs = 205 ♦ mOOO = 0.025625 f 



DTMF 
frequency 



K(floating 
point) 



K( integer) 



Absolute error 
ofK 



ofK 



697 Hz 


17.861 


18 


0.139 


0.78% 


770 Hz 


19.731 


20 


0.269 


136% 


852 Hz 


21.833 


22 


0.167 


0.76% 


941 Hz 


24.113 


24 


0.113 


0.47% 


1209 Hz 


30.981 


31 


0.019 


0.06% 


1336 Hz 


34.235 


34 


0.235 


0.69% 


1477 Hz 


37.848 


38 


0.152 


0.40% 


1633 Hz 


41.846 


42 


0.154 


037% 



The transfer function for Goertzel algorithm is 



N(1X N(9>=172; 
N(2), N(10)=177; 
N(3), N(U>=178; 
N(4XN(U>=178; 

N(5),N(13H72; 
N(6X N(14)=168 
N(7),N(13>=168 
N(8),N(16)=176 



1 



40 



l-2cos(2nJWV)r l + r 4 



where 



45 



50 



As shown in FIG. 3, the operation of the DSP is repre- 
sented as 16 Goertzel DFT blocks labeled GDFT 
(1)-GDFT — (16). Each of the Goertzel GDFT blocks, 
GDFT (1)-GDFT (16) produces a respective value A(l)-A 
(16). The values A(l)-A(16) are comprised of four smaller 55 
sub-arrays. The values A(l)-A(16) are each provided to 
respective multipliers which multiplies the respective A(n) 
value with a respective gain. These multipliers correspond to 
the gain adjustment in step 126 of FIG. 2. 

Referring now to FIG. 4, a diagram illustrating operation 
of the Goertzel algorithm is shown. As noted above, the 
Goertzel algorithm evaluates the DFT of input data with 
reduced computation and with increased efficiency. Opera- 
tion of the Goertzel algorithm is similar to a filter 
implementation, wherein the Goertzel algorithm does not 
require a buffer of input data items prior to operation, but 65 
rather computes a new output result with each occurrence of 
a new input sample. As shown in FIG. 4, the Goertzel 



Based on the above information, it is very important to 
reduce the error of the K values. 

As discussed above, the system and method of the present 
invention uses different N values for different frequencies. 
The table below describes the different N values or frame 
sizes used according to the preferred embodiment of the 
invention for the different DTMF tone frequencies. The table 
also illustrates the associated K errors: 



60 













% 


Std DTMF 




K(flcatiiig 


K(inte- 


Absolute 


error of 


frequency 


N 


point) 


ger) 


error of K 


K 


697 Hz 


172 


14585 


15 


0:014 


0.09% 


770 Hz 


177 


17.036 


17 


0D36 


0.20% 


852 Hz 


178 


18557 


19 


0.043 


0.20% 


941Hz 


178 


20537 


21 


0.063 


030% 


1209 Hz 


172 


25593 


26 


0.007 


0.02% 


1336 Hz 


168 


28.056 


28 


0.056 


0.20% 


1477 Hz 


168 


31,017 


31 


0.017 


0.05% 


1633 Hz 


176 


35526 


36 


0.074 


0.20% 



Using this method, the present invention reduces the K 
error to less than 0.3%, as compared with up to 1.36% for 
prior art methods using fixed N. In addition, the N values or 
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frame sizes are reduced from 205 to 178 or less, which The gains were computed and also verified using labora- 
allows earlier detection by 15% or more. Therefore, the DSP tory method. The gain adjustments used in the preferred 
106 preferably uses up to 178 digital speech samples which embodiment are listed below: 
are input to the Goertzel algorithm. A different number of 
these samples are used for the different tone frequencies. 5 
Application of the Goertzel algorithm produces an array of 
16 frequency domain values, wherein the outputs are basi- 
cally discrete Fourier transform results. Eight of these values 
are for the fundamental frequencies and eight are for the 
second harmonics. This is illustrated in HGS. 3 and 4. 10 

Most DFT/Goertzel based algorithms have difficulty 
meeting the required frequency distortion allowance. Hie 
frequency distortion allowance requires that all DTMF sig- 
nals with a certain percentage of frequency distortion have 
to be detected as valid DTMF tones. However, if the 15 

distortion exceeds a certain larger percentage, then detection 111 ^ P referred embodiment, the DTMF detector 102 
must be denied The frequency distortion allowance criteria does not use ^ second harmonic outputs for computing 
is described below. twist Therefore, gain adjustment for the second harmonics 

Let Fdtmf be a standard frequency and F is the actual * s * ess important. Thus, in the preferred embodiment, gain 
frequency. Then: 20 adjustment for the second harmonics is not performed in 

order to save computational resources. Assuming the gain 
if Fdtmj*(i-2S%)<F<F*mHi + 2 M), then it is a valid tone factor for frequency Fl is Gl, then the gain adjustment is 

if F>Fdtmf*(i+3.s%) or F<Fdatf*(i-3.5%) then it is not a valid calculated by performing: 



FreqfHz): 


Oaia factor 


697 


1.00000000 


770 


098902064 


852 


057439963 


941 


057849229 


1209 


1.00475907 


1336 


1.03484721 


1477 


1.02940728 


1633 


059550578 



tone. 



25 AXlXn+AQ) 



One advantage of using variable frame sizes is that the 

DTMF detector 102 more easily meets the required fre- ^ & performed on the first eight elements in the Goertzel 

quency distortion allowance. If the desired frequency is outout ^ A[1? _ , . ^ ^ &st ham ionics, and a new 

located exactly at the bin (the detecting frequency) of the array A - [lf . . . , i<q is generated. The new array comprises 

Goertzel algorithm, then +2 JS% error in the input frequency 30 me gam adjusted Goertzel DFT output It is noted that the 

decreases the output the same amount as a -2.5% error in the array elements A[9 . . . 16] are not adjusted and thus these 

signal. The reason is that both errors are at the same distance values remain the same, i.e., A*[9 . . . 16]=A[9 ... 16], 

from the bin. This is not the case if the desired frequency has step 128— Determining the Maximum Gain in Each Fre- 

error relative to the bin, it being noted that error is unavoid- 35 quency Group 

1C / 3 *?i d ^ F ° r CXample ' aSSUmC As noted above, the gain adjustment performed in step 

the desired frequency is 2.0% bigger than the frequency bin. 126 produces four sub-arrays of the first and second har- 

In this example, -2.0% error in the input signal moves the monies of the row and column frequencies of the DTMF 

frequency nearer to the bin and generates the largest outout standard frequencies. In other words, the array A*[l, , . . ,16] 

result. On the other hand, +2.5% error moves the frequency 40 is comprised of four smaller arrays. Sub-array A'[l, . . . ,4] 

further away from the bin and generates a very small output comprises the values for the fundamental row frequencies 

The asymmetric effect makes it very difficult to set the and detennines the row to which the pushed key belongs. 

threshold for frequency distortion allowance test, therefore Sub-array A'[5 8] comprises the values for the funda- 

not meeting the very stringent telecom standards. Significant mental column frequencies and determines the column to 

computational resources are necessary to determine if an 45 which the pushed key belongs. The sub-array A'[9, . . . ,12] 

energy value is at the right side or left side of the frequency comprises the DFT values for the second harmonics of the 

bin. Therefore, setting different thresholds for left and right row frequencies, and the sub-array A'[13, . . . ,16] comprises 

side distortion, i.e,, actual frequency greater or smaller than toe values for the second harmonics for the column frequen- 

the desired frequency, is not an adequate solution. 50 aes * 

Step 126 — Gain Adjustment Referring now to FIG. 5, a block diagram is shown 

In step 126 the DSP 106 preferably performs gain illustrating operation of step 128 in FIG. 2. As shown in FIG. 

adjustment, as shown in FIG. 3. Calculating the frequency 5, operation of the DSP 106 is represented as four blocks 

spectrum using different frame lengths N produces differing 151, 152, 153, and 154, wherein each of the blocks deter- 

energy contents in a plurality of the energy values. Thus the 55 niines the value in the respective sub-array that has the 

DTMF detector 102 preferably multiplies a gain value with maximum value or gain. The blocks 151 and 152 for the two 

each of the energy values to adjust the gain of each of the first harmonic sub-arrays also obtain the respective index of 

energy values. It is noted that the present invention provides me maximum gain value. The blocks 151 and 152 for the 

better results even without gain adjustment However, gain sub-arrays A' [1, ... ,4] and A'[5, ... ,8] each provide outputs 

adjustment is desirable since different frame lengths have 60 comprising the maximum value in the sub-array and the 

been applied to different frequencies and also because respective index. The blocks 153 and 154 for the sub-arrays 

different Goertzel algorithm transfer functions have different A'P* . . . ,12] and A[13, . . . ,16] provide outputs comprising 

gains at the narrow pass band. Thus the present invention toe respective maximum values in each of the sub-arrays, 

preferably adjusts the gain in order to more correctly evalu- Thus in step 128 the DSP 106 determines the maximum 

ate the results in the array A[l . . . 16]. After the gain 65 value in the first two sub-arrays and the respective indices 

adjustment in step 126, the result is a modified frequency within the sub-array. This step can be described using 

spectra values A'(m). pseudo code: 
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M[1}=0;I[1]=0; 
foi(i=l to4){ 

if A'[i>M[l] then 

M[l]=A'[i]; 

i[i] = i; 

endif 

} 

Mf2] = 0;I[2]=0; 
Cor (i = 5 to 8) { 

ifA'[i]>MI2Jthen 

Mt2]=A'[i]; 

H2] = i; 

calif 

} 



Thus, in step 128 the DSP 106 obtains the maximum 
values of the row and column frequency group (M[l], M[2]) 
and their respective index (I[l], I[2]) for the first harmonic 
sub-arrays. In the preferred embodiment, the DSP 106 does 
not determine the indices for the second harmonics. Rather 
the DSP 106 only obtains the maximum frequency output for 
the second harmonics of the row frequency, M[3J, and the 
maximum frequency output for the second harmonics of the 
column frequency, M[4], as shown in FIG. 5. 
Steps 130 and 132 — Static and Dynamic Thresholding 

The present invention also includes a method for improv- 
ing the DTMF detector's functional dynamic range and 
noise immunity by performing static thresholding followed 
by dynamic thresholding. This method increases the func- 
tional input signal dynamic range and has greater speech/ 
noise immunity, i.e., is more able to avoid detection trig- 
gered by speech or noise. The static and dynamic 
thresholding performed in steps 130 and 132 eliminate 
invalid DTMF tones based on both signal level and signal/ 
noise ratio. The use of both static and dynamic thresholding 
allows a small static threshold to be selected, thus providing 
the DTMF detector with a wider dynamic range. The 
dynamic thresholding in step 132 essentially performs a 
signal/noise ratio estimate, preferably using the ratio 
between the maximum value in the group and the second 
largest value in the sub-array. The dynamic threshold per- 
formed in step 132 provides a threshold comparison using an 
acceptable signal/noise ratio, which effectively prevents 
noise or speech-triggered detection, i.e., provides better 
speech immunity. 

Referring now to FIG. 6, a flowchart diagram illustrating 
operation of the static thresholding step DO is shown in FIG. 
1. As shown, the maximum values from the first two 
sub-arrays corresponding to the first harmonic row and 
column frequency values are received in the static thresh- 
olding step 130. As shown, in step 242 the DSP 106 
determines if M(l)>Ts, i.e., if the maximum value of the first 
sub-array is greater than Ts. If not, then no detection occurs 
and operation completes. If M(l)>Ts in step 242, then in step 
244 the DSP determines if the value M(2) is greater than Ts. 
If not, then again no detection is determined. If both M(l) 
and M(2) are greater than the threshold value Ts, then 
dynamic thresholding is performed on these values in step 
132 of FIG. 2. 

Referring now to FIG. 7, the dynamic thresholding step 
132 is shown. As shown, the dynamic thresholding step 
compares the ratio of the maximum value M(l) in the first 
sub-array to each of the other values in the sub-array A'[l], 
A'[2], A'[3] and A'[4], and sets the signal/noise ratio error as 
the minimum ratio. Likewise, the DSP 106 in step 132 
determines the minimum ratio of M(2) with each of the other 
values from the sub-array A'[S], A [6], A'[7] and A'[8], and 
sets the signal/noise ratio accordingly. More specifically, in 
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FIG. 7 the DSP 106 compares the ratio of the maximum 
value Ml and each of the other values in the sub-array A'[l], 
A'[2] , A'[3] and A'[4] with a value Td and sets the SNR error 
as the ratio which is less than Td. likewise, the DSP 106 

5 compares the ratio of the maximum value M(2) and each of 
the other values in the sub-array A'[5], A'[6], A'[7] and A*[8] 
with a value Td and sets the SNR error as the ratio which is 
less than Td Because Td=10 Cns9Uim/ SN * * d * /20) , setting the 
SNR error means that either or both the column and row 

10 frequency values do not meet the required SNR. If SNR 
error is set, no detection is indicated. 

Thus the frequency domain values (A'[l, . . . ,16]) 
obtained through the Goertzel algorithm and gain adjust- 
ment are compared against thresholds using both static and 

15 dynamic techniques. In the preferred embodiment, the static 
threshold is set to 6 (which provides a theoretical functional 
dynamic range of over 50 dB) and the dynamic threshold is 
set to 4 (about 12 dB). This allows very small DTMF signals 
(wide dynamic range) to be detected, as long as the signals 

20 have an acceptable signal to noise ratio. This also effectively 
eliminates most idle channel noise, speech, and white noise. 
Step 134 — Second Harmonic Thresholding 

In step 134, the DSP 106 performs second harmonic 
thresholding to further examine the validity of the detected 

25 tones and eliminate speech triggered DTMF detection. This 
occurs before DTMF detection and usually is a continuous 
task for the DSP 106, even after the channel has been 
established The human voice, especially the female voice, 
can have a large amount of energy components over 1000 

30 Hz and can be mistaken as DTMF signal. Some music can 
also trigger DTMF detection. Thus speech detection and/or 
speech immunity is a very important criteria for evaluating 
the quality of a DTMF detector. 
As discussed above, a real DTMF signal is the sum of two 

35 sinusoids and has two steep peaks in the frequency domain. 
Thus, a real DTMF signal does not have significant energy 
at second or higher harmonics. Speech on the other hand 
generally always has a significant amount of energy at the 
second and higher harmonics. This characteristic of speech 

40 make it easier to distinguish speech from DTMF signal. 
Thus the present invention examines the second harmonics 
of the fundamental DTMF frequency. 

Prior art methods have conventionally examined the value 
of the second harmonics, A'[9, . . . ,16] ■ According to the 

45 present invention, the DSP 106 compares the second har- 
monics with the values at the fundamental harmonics, A'[l, 
. . . ,8]. The DSP 106 computes the ratio of M[l]/M[3] and 
M[2]/M[4], and compares the results with a pre-set thresh- 
old 

50 It is noted that the ratio A'[i]/A'[i+8] (assuming A'[i] is the 
biggest value in its group) is not used because, if the speech 
has some energy at a frequency x% away from the bin, its 
second harmonic will have a distance of 2x% from the 
second harmonic bin. Therefore, the present invention uses 

55 the maximum value in the second harmonic group to provide 
an estimate of the energy at the second harmonics. This is 
beneficial to talk-off performance of the DTMF detector. 

Referring now to HG. 8, a flowchart diagram illustrating 
operation of the second harmonic thresholding performed in 

60 step 134 of FIG. 2 is shown. As shown, in FIG. 8 the DSP 
106 compares the ratio of the maximum values of the first 
and second harmonics, M1/M3, with the threshold value 
Tsh, and also compares the ratio of the maximum values of 
the first and second harmonics of the column frequency, Le., 

65 the ratio M2/M4, and determines if this value is greater than 
the threshold Tsh. If either of the ratios M1/M3 or M2/M4 
are not greater than the threshold Tsh, then a speech error 
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signal is generated indicating that the detected signal actu- still passes the guard-time check. If this comparison fails 
ally corresponds to speech and not to a DTMF signal. again, then the signal on the current frame fails the guard- 
Step 136— -Guard Time Check time check and cannot be recognized as a valid DTMF 

The present invention discloses an improved guard time signal, 

check method which overcomes the difficulty of meeting the 5 step 13 8— Twist Computation and Thresholding 

guard-time performance requirement when using frame jf the signal passes all of the above steps, then the final 

based processing methods (such as the Goertzel algorithm). stage of eva i aat i on i s me twist check. Like the other param- 

The present invention also uses a larger cut-off duration eters described above, twist is another important criteria to 
value (but less than Tl) as described below. A larger cut-off eliminate unqualified signals. The twist of a DTMF signal is 

duration value provides better talk-off performance for the 10 defined as the level difference in dB between the higher and 

DTMF detector. This is because speech triggered fake ^ lower tone. 
DTMF signals' tend to have a short duration, while also 

having a frequency domain energy presentation which ^(2) 

appears like a DTMF signal. The guard-time parameter Twtst=20iog I0 ^ 

requires that all DTMF signals having a duration longer than l s 

Tl be detected as valid DTMF tones, i.e., no misses are . . , . . A 

..... .. , I . ™ Thus, if the twist is positive, the higher frequency level tone 

allowed), whereas all DTMF signals shorter than T2 /4 , y . . r . . ' fr *u i c 

rr>xri w*n ™* a^^a « iytmij c *„«*i« ah ( me column frequency) is greater than the lower frequency 

(ira) can not be detected as valid DTMF ignak. AU 1 d tone ^ ^ frfi } m . Q ^ > 

tones have a duration between .Tl an T2 axe ^n't care ^ ff \ t J" is neg ative, i.e., if the higher 

which means the detector has the flexibility to either detect 20 , . . , , ° \ . \ ^ 

. * u . . T . * j j * 77 ^7 frequency level tone (the column frequency) is less than the 

or miss the signal. In general, it is deemed desirable to have , ~ , \ „ /4 . N ^ 

\*u „i *: -c lower frequency level tone (the row frequency), then the 

a fixed cut-off value to warrant the guard-time performance. ^ SAJ « JL , . ^ » « \ j ^ . ,, 

t*u at clt * j a a * *i * / . . twist is called backward twist A telecom standard typically 

The AT&T standard, and most other standards, require . ~ *. • * n u „ * ~i *. ■ * *iT u ,i 

. ryrijiT? t *u ,in * u defines a twist allowance between a forward twist threshold 

that any DTMF sequence longer than 40 ms must be . , . . u , . „ , 

~ * a ** * ii *u v • a ivpiic « and a backward twist threshold. For example, AT&T 

accepted if it meets all other criteria. A DTMF sequence 25 - eflU « res . 

shorter than 20 ms should never be detected. Since the ^ 

DTMF detector 102 of the preferred embodiment processes _g dB<twist< ^ 4 & 

the speech sample as blocks, the DTMF detector 102 com- K 

ii^.i . Af\ channel, distortion and noise are added to the DTMF signal 

signal lasts at least 40 ms. . , , . . . , to .. 

„ r . . n a . , .„ ^ . and cause phase delay, group delay and even severe ampli- 

Refemng now to FIG. 9, a flowchart diagram illustrating . . .. , *\ r~ Jr * v * ^ - . ^ 

. 2_ . , ' . , j • * * tude distortion. These factors complicate the twist compu- 

the guard time check operation performed in step 136 of . A * • *u ivmirj jT 

HG 1 is shown. As shown, in step 312 the DSP 106 ^ ^J** to *f f™ 6111 ™ on t ™ d ? B0 : 

dctennines if the index values I[l] and are equal and 35 ' OT 102 **** ™ f*??** 11 ^ w ^ e » * e r ^ nl 4 

index values I[2] and T[2] are equal. Ttai. this ste^, deter- S^Tt^I^^^^T^JV^ « t ? a "J 

mkesiftheimiFsignalhaslastedac^rtatoperiodof^ M[2], which are the gain adjusted outputs from step 128, and 

More specifically, thf step determines if Eta£2rf only update the lvalue when both the^er and lower 

maximum values in consecutive frames are identical. If so, " ne level m * e frame » ^ * e ? ^ P"! 10 " 5 

~TT . „t " . . i frame, assumine that the row and column indexes. Iffl and 

then in step 314 the DSP 106 computes the ratio or gain of 40 t£i ' ™». ,u ™ , » UMl "^r . w rrrr" 7.. ,. ™ , J: 

me maximum value from the c urrent and prior frame for »™in Jta same. This method is dlustated m FIG. 10 

each of the row and column frequencies. In other words, in aadu executed 88 ^ foUowm 8 P^do-code: 

step 314 the DSP 106 computes: 



nn-wfivMYiv m] = r[i] AND 112] = T[2] THEN 

m»=«(l/M(i), 45 IFM[l>M'Il]ANDM[2)>M't2]THHN 

0(2^(2X^X2); ^ ^ = 20"o g 10(M[ 2V M[l]) 

TWist — TWist* 

In steps 316 and 318 the DSP 106 delOTnines if the e^otf 

computed ratios Gl and G2 are within a range >1/TG and endif 

<TG. so — 

If the DSP 106 detects a DTMF signal which meets all of T . t . , . f , . 

the criteria described above in steps 130-134, then the DSP * F"? cal "nplfmentation, only a linear twist value is 

106 compares the maximum values with the values at the computed to avoid logarithmic computation, and the linear 

previous frame only if the indices remain the same. If the twist value IS ^ com P ared Wlth ^ thresholds to 

ratio, i.e. f the result of the division, of the rnaximum values 55 4**™*™ tf me detection meets the twist requirement Note: 

fromthecurrentandr^iorframeis greater than the value Tg, Twtst 18 me ^ vaiue for ^ Prev 1008 
then the signal passes the guard-time check. It is noted that 

the value Tg is a programmable value that can be changed T[1 ]^ ] ' m ^ ^ r[2] * m m ** pmvkms 
to meet different guard-time requirements. The value Tg is 

^erably 2^ Thus in ^ step 136 of FIG 2 the DTMF 60 M[1] ^U*** p*^ ^ M[2] is M[2i in the pre- 

detector 102 determines if the detected signal is of sufficient vious frame 
length to actually be a DTMF signal. 

If the guard time check fails, the DSP 106 waits for the DTMF Receiver Results Using MITEL 7292 Test 

computation of next frame and compares the maximum Tape 
value of current frame with the maximum value in the next 65 

frame, also assuming that the row and column indices The following tests were performed on the DTMF detec- 

remain the same. If the result is greater than Tg, the signal tor of the present invention using the Mitel 7292 test tape. 
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Tfcst 1: Frequency Distortion Allowance 
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Test 3: Dynamic Range Test 

34- dB (Maximum Input Signal dynamic range) 
Test 4: Guard-time Test: 28.9 ms. 
Test 5: Signal to Noise Ratio Test 

Input SNR 24 dB 18 dB 12 dB 
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Test 6: Talk-off Test 

6 hits (one real DTMF signal and rive raise detections). 
System Requirements: 

Processing power <1.2 MIPS. 

ROM usage; 314 words. 

RAM usage: 69 words. 



Conclusion 

Therefore, the present invention comprises a system and 
method fox performing DTMF detection with improved 
speed and accuracy. The combined methods of using varying 
frame sizes, combined static/dynamic thresholding, second 
harmonic thresholding, guard-rime protection, and twist 
computation provide very good talk-off performance. The 
DTMF detector 102 of the present invention only has about 
6 hits running the Mitel CM7291 standard test tape, whereas 
an algorithm which registers less than 10 hits is considered 
very good. The present invention also provides excellent 
performance in meeting other parameter criterias. 

Although the system and method of the present invention 
has been described in connection with the preferred 
embodiment, it is not intended to be limited to the specific 
form set forth herein, but on the contrary, it is intended to 
cover such alternatives, modifications, and equivalents, as 45 
can be reasonably included within the spirit and scope of the 
invention as defined by the appended claims. 

We claim: 

1. A method for detecting dual tone multifrequency 
(DTMF) signals, comprising: 
receiving a plurality of digital samples of a received 
signal, wherein said received signal includes a plurality 
of tones, wherein said plurality of tones comprise two 
or more tones from a plurality of different uncorrelated 
frequencies, wherein said plurality of different uncor- 
related frequencies comprise two or more frequency 
groups; 

calculating a frequency spectrum on said plurality of 
digital samples for each of said plurality of different 
uncorrelated frequencies, wherein said calculating uses 60 
a frame length N comprising at least a subset N of said 
digital samples for each of said different uncorrelated 
frequencies, wherein said calculating uses a different 
frame length N of said digital samples for at least a 
plurality of said different uncorrelated frequencies, 
wherein said calculating produces an energy value for 
each of said different uncorrelated frequencies; 



50 



55 



determining maximum values of said energy values for 
each of said two or more frequency groups to detect 
said plurality of tones in said received signaL 

2. The method of claim 1, wherein said calculating a 
frequency spectrum using different frame lengths N pro- 
duces differing energy contents in a plurality of said energy 
values, the method further comprising: 

multiplying gain values with at least a plurality of said 
energy values to adjust the gain of said energy values, 
wherein said multiplying occurs after said calculating a 
frequency spectrum on said plurality of digital samples 
and prior to said determining maximum values of said 
energy values. 

3. The method of claim 2, wherein said multiplying gain 
values with at least a plurality of said energy values to adjust 
the gain of said energy values comprises multiplying dif- 
ferent gain values with at least a plurality of said energy 
values, wherein the magnitudes of said different gain values 
depend on the size of said different frame lengths N for said 
different uncorrelated frequencies, wherein the magnitudes 
of said different gain values further depend on the difference 
between said uncorrelated frequencies and said calculated 
frequency spectrum; 

wherein said different gain values compensate for said 
differing energy contents in said plurality of said energy 
values. 

4. The method of claim 1, wherein said different frame 
lengths N are designed to optimally align said calculated 
frequency spectrum at said different uncorrelated frequen- 
cies for each of said plurality of different uncorrelated 
frequencies. 

5. The method of claim 1, wherein said calculating a 
frequency spectrum on said plurality of digital samples for 
each of said plurality of different uncorrelated frequencies 
comprises performing a Goertzel algorithm on said plurality 
of digital samples for each of said plurality of different 
uncorrelated frequencies; 

wherein said performing said Goertzel algorithm uses a 
different frame length N of said digital samples for a 
plurality of said different uncorrelated frequencies. 

6. The method of claim 1, wherein said calculating a 
frequency spectrum on said plurality of digital samples for 
each of said plurality of different uncorrelated frequencies 
comprises performing a Discrete Fourier Transform (DFT) 
on said plurality of digital samples for each of said plurality 
of different uncorrelated frequencies; 

wherein said performing said DFT uses a different frame 
length N of said digital samples for a plurality of said 
different uncorrelated frequencies. 

7. The method of claim 1, wherein said plurality of 
different uncorrelated frequencies comprise two frequency 
groups including a first frequency group comprising fre- 
quencies 697 Hz, 770 Hz, 852 Hz, and 941 Hz, and a second 
frequency group comprising frequencies 1209 Hz, 1336 Hz, 
1477 Hz, and 1633 Hz. 

8. The method of claim 7, wherein said different frame 
lengths N have the following values for said first and second 
frequency groups: 

DTMF frequency N 



65 



697 Hz 
770 Hz 
852 Hz 
941 Hz 



172 
177 
178 
178 
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-continued 



DTMF frequency 


N 


1209 Hz 


172 


1336 Hz 


168 


1477 Hz 


168 


1633 Hz 


176. 



9. The method of claim 1, wherein said calculating a 
frequency spectrum comprises calculating a frequency spec- 
trum on said plurality of digital samples for first and second 
harmonics of each of said plurality of different uncorrelated 
frequencies; 

wherein said calculating uses the same frame length N of 
said digital samples for both said first and second 
harmonics, and uses a different frame length N for a 
plurality of said different uncorrelated frequencies. 

10. The method of claim 1, further comprising: 
receiving said received signal, wherein said received 

signal includes a plurality of tones, wherein said plu- 
rality of tones comprise two or more tones from a 
plurality of different uncorrelated frequencies, wherein 
said plurality of different uncorrelated frequencies 
comprise two ox more frequency groups; and 
performing analog to digital conversion on said received 
signal to produce said plurality of digital samples of 
said received signal. 

11. A dual tone multifrequency detector for detecting dual 
tone multifrequency (DTMF) signals, comprising: 

means for receiving a plurality of digital samples of a 
received signal, wherein said received signal includes a 
plurality of tones, wherein said plurality of tones com- 
prise two or more tones from a plurality of different 
uncorrelated frequencies, wherein said plurality of dif- 
ferent uncorrelated frequencies comprise two or more 
frequency groups; 

a digital signal processor for calculating a frequency 
spectrum on said plurality of digital samples for each of 
said plurality of different uncorrelated frequencies, 
wherein said digital signal processor uses a frame 
length N comprising at least a subset N of said digital 
samples for each of said different uncorrelated 
frequencies, wherein said digital signal processor uses 
a different frame length N of said digital samples for at 
least a plurality of said different uncorrelated 
frequencies, wherein said digital signal processor pro- 
duces an energy value for each of said different uncor- 
related frequencies; 

wherein said digital signal processor determines maxi- 
mum values of said energy values for each of said 
two or more frequency groups to detect said plurality 
of tones in said received signal. 

12. The dual tone multifrequency detector of claim 11, 
wherein said digital signal processor calculating a frequency 
spectrum using different frame lengths N produces differing 
energy contents in a plurality of said energy values; 

wherein said digital signal processor multiplies gain val- 
ues with at least a plurality of said energy values to 
adjust the gain of said energy values. 

13. The dual tone multifrequency detector of claim 12, 
wherein said digital signal processor multiplies different 
gain values with at least a plurality of said energy values, 
wherein the magnitudes of said different gain values depend 
on the size of said different frame lengths N for said different 
uncorrelated frequencies; 

wherein the magnitudes of said different gain values 
further depend on the difference between said uncor- 
related frequencies and said calculated frequency spec- 
trum; 
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wherein said different gain values compensate for said 
differing energy contents in said plurality of said energy 
values. 

14. The dual tone multifrequency detector of claim 12, 
5 wherein said different frame lengths N are designed to 

optimally align said calculated frequency spectrum at said 
different uncorrelated frequencies for each of said plurality 
of different uncorrelated frequencies. 

15. The dual tone multifrequency detector of claim 11, 
10 wherein said digital signal processor performs a Goertzel 

algorithm on said plurality of digital samples for each of said 
plurality of different uncorrelated frequencies; 
wherein said digital signal processor uses a different 
frame length N of said digital samples for a plurality of 
15 said different uncorrelated frequencies in performing 
said Goertzel algorithm. 

16. The dual tone multifrequency detector of claim 11, 
wherein said digital signal processor performs a Discrete 
Fourier Transform (DFT) on said plurality of digital samples 

20 for each of said plurality of different uncorrelated frequen- 
cies; 

wherein said digital signal processor uses a different 
frame length N of said digital samples for a plurality of 
said different uncorrelated frequencies in performing 
25 said DFT. 

17. The dual tone multifrequency detector of claim 11, 
wherein said plurality of different uncorrelated frequencies 
comprise two frequency groups including a first frequency 
group comprising frequencies 697 Hz, 770 Hz, 852 Hz, and 
941 Hz, and a second frequency group comprising frequen- 
cies 1209 Hz, 1336 Hz, 1477 Hz. and 1633 Hz. 

18. The dual tone multifrequency detector of claim 17, 
wherein said different frame lengths N have the following 
values for said first and second frequency groups: 
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697 Hz 


172 


770 Hz 
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19. The dual tone multifrequency detector of claim 11, 
wherein said digital signal processor calculating a frequency 
spectrum comprises calculating a frequency spectrum on 

50 said plurality of digital samples for first and second har- 
monics of each of said plurality of different uncorrelated 
frequencies; 

wherein said calculating uses the same frame length N of 
said digital samples for both said first and second 
55 harmonics, and uses a different frame length N for a 
plurality of said different uncorrelated frequencies. 

20. A dual tone multifrequency detector for detecting dual 
tone multifrequency (DTMF) signals, comprising: 

a coder/decoder (codec) which receives an analog input 
60 signal and generates a plurality of digital samples of the 
received analog input signal, wherein said received 
analog input signal includes a plurality of tones, 
wherein said plurality of tones comprise two or more 
tones from a plurality of different uncorrelated 
65 frequencies, wherein said plurality of different uncor- 
related frequencies comprise two or more frequency 
groups; 
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a digital signal processor for calculating a frequency 
spectrum on said plurality of digital samples for each of 
said plurality of different uncorrelated frequencies, 
wherein said digital signal processor uses a frame 
length N comprising at least a subset N of said digital 
samples for each of said different uncorrelated 
frequencies, wherein said digital signal processor uses 
a different frame length N of said digital samples for at 
least a plurality of said different uncorrelated 
frequencies, wherein said digital signal processor pro- 
duces an energy value for each of said different uncor- 
related frequencies; 

wherein said digital signal processor determines maxi- 
mum values of said energy values for each of said 
two or more frequency groups to detect said plurality 
of tones in said received signal. 

21. A method for detecting multi tone nuiltifrequency 
(MTMF) signals, comprising: 

receiving a plurality of digital samples of a received 
signal, wherein said received signal includes a plurality 
of tones, wherein said plurality of tones comprise two 
or more tones from a plurality of different uncorrelated 
frequencies, wherein said plurality of different uncor- 
related frequencies comprise one or more frequency 
groups; 

calculating a frequency spectrum on said plurality of 
digital samples for each of said plurality of different 
uncorrelated frequencies, wherein said calculating uses 
a frame length N comprising at least a subset N of said 
digital samples for each of said different uncorrelated 
frequencies, wherein said calculating uses a different 
frame length N of said digital samples for at least a 
plurality of said different uncorrelated frequencies, 
wherein said calculating produces an energy value for 
each of said different uncorrelated frequencies; 

determining maximum values of said energy values for 
each of said one or mare frequency groups to detect 
said plurality of tones in said received signal. 

22. Hie method of claim 21, wherein said calculating a 
frequency spectrum using different frame lengths N pro- 
duces differing energy contents in a plurality of said energy 
values, the method further comprising: 

multiplying different gain values with at least a plurality 
of said energy values to adjust the gain of said energy 45 
values, wherein said multiplying occurs after said cal- 
culating a frequency spectrum on said plurality of 
digital samples and prior to said determining maximum 
values of said energy values. 

23. The method of claim 21, wherein said different frame 
lengths N are designed to optimally align said calculated 
frequency spectrum at said different uncorrelated frequen- 
cies for each of said plurality of different uncorrelated 
frequencies. 

24. The method of claim 21, wherein said calculating a 
frequency spectrum on said plurality of digital samples for 
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each of said plurality of different uncorrelated frequencies 
comprises performing a Goertzel algorithm on said plurality 
of digital samples for each of said plurality of different 
uncorrelated frequencies; 
wherein said performing said Goertzel algorithm uses a 
different frame length N of said digital samples for a 
plurality of said different uncorrelated frequencies. 

25. A multi tone multifrequency detector for detecting 
multi tone multifrequency (MTMF) signals, comprising: 

means for receiving a plurality of digital samples of a 
received signal, wherein said received signal includes a 
plurality of tones, wherein said plurality of tones com- 
prise two or more tones from a plurality of different 
uncorrelated frequencies, wherein said plurality of dif- 
ferent uncorrelated frequencies comprise one or more 
frequency groups; 

a digital signal processor for calculating a frequency 
spectrum on said plurality of digital samples for each of 
said plurality of different uncorrelated frequencies, 
wherein said digital signal processor uses a frame 
length N comprising at least a subset N of said digital 
samples for each of said different uncorrelated 
frequencies, wherein said digital signal processor uses 
a different frame length N of said digital samples for at 
least a plurality of said different uncorrelated 
frequencies, wherein said digital signal processor pro- 
duces an energy value for each of said different uncor- 
related frequencies; 

wherein said digital signal processor determines maxi- 
mum values of said energy values for each of said 
one or more frequency groups to detect said plurality 
of tones in said received signal. 

26. The multi tone multifrequency detector of claim 25, 
wherein said calculating a frequency spectrum using differ- 
ent frame lengths N produces differing energy contents in a 
plurality of said energy values; 

wherein said digital signal processor multiplies a gain 
value with at least a plurality of said energy values to 
adjusfthe gain of said energy values. 

27. The multi tone multifrequency detector of claim 25, 
wherein said different frame lengths N are designed to 
optimally align said calculated frequency spectrum at said 
different uncorrelated frequencies for each of said plurality 
of different uncorrelated frequencies. 

28. The multi tone multifrequency detector of claim 25, 
wherein said digital signal processor performs a Goertzel 
algorithm on said plurality of digital samples for each of said 
plurality of different uncorrelated frequencies; 

wherein said digital signal processor uses a different 
frame length N of said digital samples for a plurality of 
said different uncorrelated frequencies in performing 
said Goertzel algorithm. 
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